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Abstract 

We present the quantum phase diagram of mesoscopic SQUID array. We predict different quan- 
tum phases and the phase boundaries along with several special points. We present the results of 
magnetic flux induced dissipation on these points and also on the phase boundaries. The Josephson 
couplings at these points are dependant on the system parameter in a more complicated fashion 
and differ from the Ambegaokar and Baratoff relation. We derive the analytical relation between 
the dissipation strength and Luttinger liquid parameter of the system. We observe some interesting 
behaviour at the half-integer magnetic flux quantum and also the importance of co-tunneling effect. 
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1. INTRODUCTION 



Josephson junction arrays have attracted considerable interest in the recent years owing 
to their interesting physical properties like quantum phase transition, quantum critical be- 
haviour and Coulomb blockade etc The most universally observed behaviour for 
the low- dimensional superconducting system is the superconductor-insulator transition in 
superconducting Aim, wires, Josephson junction array in one and two dimensional giving 
rise to intense experimental and theoretical activity 3Hl7| This superconductor-insulator 
transition occurs at low temperature ( mili-Kelvin scale) due the variation of one of the pa- 
rameter of the system such as the normal state resistance of the Aim, thickness of the wire, 
the Josephson coupling of Josephson junction array and the magnetic flux of the mesoscopic 
SQUID array 

In one dimensional superconducting quantum dot lattice system (mesoscopic SQUID array, 
Cooperpair box array) one elementary excitation occurs, i.e., the quantum phase slip cen- 
ters (QPS) lo|, 14, Q, Q- Recently the physics of QPS and the related physics has been 
studied extensively. Quantum phase slip process is a topological excitation, it's a discrete 
process in space-time domain of a one dimensional superconducting system. In this process 
the amplitude of the order parameter destroys temporarily at a particular point that leads 
the phase of the order parameters to change abruptly (in units of 2 tc). This process occurs 
at the time of macroscopic quantum tunneling of the system. As we understand from our 
previous study [llj and also from the existing literature that the QPS process initiate the 
superconductor- insulator transition [l^- It is well known after the work of Calderia 



and Leggett 



191 ] the dissipation plays a central role in the macroscopic quantum system. 



There are several study in the literature based on the Calderia-Leggett formalism to explain 



the different physical properties of low- dimensional tunneling j unctio n sy stem 2014261] . We 



are mainly motivated from the experimental findings of Chow et al [13j. They have done 
some pioneering work to study the effect of magnetic flux on mesoscopic SQUIDs array. 
They have observed the magnetic flux induced dissipation driven superconductor-insulator 
quantum phase transition and the evidence of qunatum critical point. The other interesting 
part of this study the length scale dependent superconductor-insulator transition and the 
magnertic flux induced superconductivity. 

The goal of this paper is to provide a detailed analysis of the appearance of different quan- 
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FIG. 1: A. Schematic diagram of one dimensional array of small capacitance dc SQUIDs. Each 
plaque is a SQUID with two Josephson junctions marked by the cross. Circle with plus sign 
represent the applied magnetic flux $. B. Equivalent representation of system A, where the dots 
are connected through tunnel junctions and the Josephson couplings of this system is tunable due 
to presence of magnetic flux <£. 

turn phases and phase boundaries in the presence and absence of dissipation for a mesoscopic 
SQUID array at T = 0. We assume the presence of ohomic dissipation motivated by the 



experiment in Ref. (|13| ) and we study the consequence of it in the different quantum 
phases and phase boundaries and also on the special points like the particle-hole symmetric 
point, charge degeneracy point, and multicritical points. We also derive the relation between 
the dissipation strength and Luttinger liquid parameter of the system. We are able to find 
out the relation between the Josephson coupling and the other interaction parameters of 
the system through the analysis of Luttinger liquid parameter of the system , we will notice 
;hat this analytical relation is more complicated and differs from the Ambegaokar-Baratoff 



271 ] relation. To the best of our knowledge, this is the first derivation in the literature. The 
plan of the manuscript is as follows. In section 2A, we introduce the basic concept of the 
appearence of quantum phase slip centers and the source of dissipation for one dimensional 
superconducting quantum dot lattice. In section 2B, we derive the relation between the 
dissipation strength and the Luttinger liquid parameter of the system. In section 3, We 
present the analysis of quantum phases and phase boundaries, in the presence and absence 
of magnetic flux induced dissipation. Section IV, is devoted for summary and conclusions. 



2. BASIC ASPECTS OF QUANTUM PHASE SLIP CENTERS AND ANALYTICAL 
DERIVATIONS OF DISSIPATIVE STRENGTHS IN TERMS OF LUTTINGER 
LIQUID PARAMETER 
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2A. Basic Aspects of Quantum Phase Slip Centers: 

It is well known to us that the Josephson coupling of the SQUID is modulated by a facror, 
Eji = Ej\cos(^)\, where 4> is the magnetic flux and </> = |^ is the magnetic flux quantum. 
Therefore we can consider the the mesoscopic SQUIDs array as a superconducting quantum 
dot (SQD) lattice with modulated Josephson junction 0, [lo|. Fig. 1 shows the equivalance 
between the mesoscopic SQUIDs array and SQD lattice with modulated Josephson junctions. 
Here we prove the appearance of QPS in SQD array with modulated Josephson coupling 
through an analysis of a minimal model. We consider two SQDs are separated by a Josephson 
junction. These two SQDs are any arbitrary SQDs of the array. Appearance of QPS is an 
intrinsic phenomenon (at any junction at any instant) of the system. The Hamiltonian of 
the system is 

H = £ ^ - Ej\cos(^-)\Y:cos(e i+1 - 6,) (1) 

where rij and 9^ are respectively the Cooper pair density and the superconducting phase of 
the i'th dot and C is the capacitance of the junction. The first term of the Hamiltonian 
present the Coulomb charging energies between the dots and the second term is nothing but 
the Josephson phase only term with modulated coupling, due to the presence of magnetic 
flux. We will see that this model is sufficient to capture the appearance of QPS in SQD 
systems. Hence this Hamiltonian has sufficient merit to capture the low temperature dissi- 
pation physics of SQD array. In the continuum limit, the partition function of the system 
is given by, Z = J D0(x,t) e~ S ^ T \ where S Q = J dr J dx^ [(«9 T #(x, r)) 2 + (d x 9(x,T)) 2 }. 
This action is quadratic in scalar-field 9(x,t), where 9(x, r) is a steady and differentiable 
field, so one may think that no phase transition can occur for this case. This situation 
changes drastically in the presence of topological excitations for which 8(x, r) is singular 
at the center of the topological excitations. So for this type of system, we express the 
9(x, t) into two components: 9(x, r) = 9q(x, t) + 9±(x, r), where 9q(x, t) is the contribution 
from attractive interaction of the system and 9i(x, r) is the singular part from topological 
excitations. We consider at any arbitrary time r, a topological excitation with center 
at X(t) = (xq(t), To(t)) . The angle is measured from the center of topological 
excitations between the spatial coordinate and the x-axis 6 l 1 (x,r) = tan~ 1 ( -^z^) ■ The 
derivative of the angle is V x 9\(x — X{r)) = n^ZxTPj^ K r o — r )' ( x ~ ^o)] which has 
a singularity at the center of the topological excitation. Finally we get an interesting 



4 



result when we integrate along an arbitrary curve encircling the topological excitations. 
J C1 dx V x 8i (x — X(t)) = 2tt. So we conclude from our analysis that when a topological 
excitation is present in the SQD array, the phase difference ,9, across the junction of 
quantum dots jumps by an integer multiple of 2ir. This topological excitation is nothing but 
the QPS in the (x, r) plane. According to the phase voltage Josephson relation ,Vj = 
(t is the time), a voltage drop occurs during this phase slip, which is the source of dissipation. 



2B. Analytical Derivations of Dissipative Strengths in Terms of Luttinger Liquid 
Parameters 

Here we derive analytical expression of magnetic flux induced dissipative strength (a) in 
terms of the interactions of the system. At first we derive the dissipative action/partition 
function of a quantum impurity system. In our study, we consider the backscattering pro- 
cess from a quantum impurity to compare the effective action with the superconducting 
tunnel junction system. We will see that the analytical structure of this dissipative action 
is identical with the mesoscopic SQUIDs array. In a different context, the authors of Ref. 



281 ] have predicted the dissipation driven quantum phase transition to occur. They have 
also considered the presence of backscattering events originated in the LL under the effect 
of dynamically screened Coulomb interaction. 

Here we consider that the impurity is present at the origin where the fermions scatter from 
the left to the right and vice versa. The Hamiltonian describing this process is 

E x = V (R\0)L(0) + h.c) =V Q [ dxd(x)cos9(x). 

The total Hamiltonian of the system is 

H = H + V [ dxd(x)cos9(x 1 r) (2) 

H = i-| uK{d x 9(x,r)) 2 + -|(^(x,r)) 2 , (3) 
The corresponding Lagrangian of the system is 



-(d T <j)(x,T)Y + u(d(f>(x,T)Y 



2ttK J u 

+V J dxS(x)cos(9(x, r)) = L + L x (4) 



where L and L\ are the non-interacting and the interacting part of the Lagrangian and K is 
the Luttinger liquid parameter of the system. The only non-linear term in this Lagrangian 
is expressed by the field 9(x = 0). We would like to express the action of the system as an 
effective action by integrating the field 9(x 7^ 0). Therefore one may consider 9(x ^ 0) as a 
heat bath, which yields the source of dissipation in the system. The constraint condition for 
the integration is 9{r) = 6{x = 0,r). We can write the partition function of the quantum 
impurity system. Z = J D9(x,t)D9(t)5(9(t) - 9(0, r))e~Xf Ldr Now we use the standard 
trick of introducing the Lagrange multiplier with auxiliary field A(r). 

Z = J D0(r)e- f° LldT J D\{T)e- iX{ - T)e( - T) 

J 'D0(x,T)ef°(- Lo+iX ™°' T » dT (5) 

The Fourier transform of the first term of Eq.(3) is 
Lq = J2 q J2iui n ^y^- 3 - ^) iw n )0(—q, —iu n ). At first we would like to calculate 
the integral: Jq dr[L — i\(r)9(0, r)], we can write this term as 

S S n 2nKv ^ iUn ^~ q > ~ iuJn "> 

(\(iu n )9(-q, -iu n ) + \(-iu n )9(q,iu n ). (6) 



2V 7 ! 

This integral appears in the integral 9(x, r). This integral is quadratic in 9. Now we would 
like to perform the Gaussian integration by completing the square. We can write the result 
as irX^ n — 2* Kv 2 2 - In the infinite length limit one can write, 

1 ttKv r dq_ ttKv _ nK 

2L ^1 ui n 2 + v 2 q 2 ~ J 2vr ui n 2 + v 2 q 2 ~ 4ui n " 

Now we would like to append this result of integration in the second integral of Z, i.e., the 
integral over A. One can write the integrand as 

S(~4 — K iuj n)K- iuj n) + -(X(iuj n )9(-q, -iu n ) 

iu>„ n 

+\(-iuj n )9(q,iuj n ). (7) 

This integral is again the quadratic integral of A, therefore, the Gaussian integral can 
be performed by completing the square. We perform the integration and finally obtain 
y, iuJn f^9(iu n )9(—icu n ). From this analytical expression, we obtain the effect of bath on 
9(t). Appearance of the factor u n signifies the dissipation. Therefore, the effective action 
reduces to 

S = £^(^n)0H<> + JdrV cos9(r) (8) 

U>Jn 



The above action implies that a single particle moves in the potential Vqcos9{t) subject to 
dissipation with friction constant , 

Now we calculate the dissipative action of mesoscopic SQUIDs array. Here, we calculate 
the effective partition function of our system. Our starting point is the Calderia-Legget fl9 ] 
formalism. Following this reference we write the action as 

Sq is the action for non-dissipative part and Si is the standard action for the system with 
tilted wash-board potential |20|, |2l| to describe the dissipative physics for low dimensional 



superconducting tunnel junctions, a is the dissipative strength of the system, 

R Q i^i / 1n x 
a = — — cos| — — | (10) 

R s O 

(the extra cosine factor which we consider in a is entirely new which probes the effect of 
an external magnetic flux and is also consistent physically), u m = is the Matsubara 

frequency and Rq (= 6A5kQ) is the quantum resistance, R s is the tunnel junction resis- 
tance, P is the inverse temperature. In one of our previous work, we have explained few 
experimental findings by using the above expression of a. In the strong potential, tunneling 
between the minima of the potential is very small. In the imaginary time path integral for- 
malism, tunneling effect in the strong coupling limit can be described in terms of instanton 
physics. In this formalism, it is convenient to characterize the profile of 9 in terms of its 
time derivative, 

^ = E^(r-r 4 ), (11) 

where h(r — Tj) is the time derivative at time r of one instanton configuration, Tj is the 
location of the i-th instanton, = 1 and —1 is the topological charge of instanton and anti- 
instanton respectively. Integrating the function over h from — oo to oo, we get drh{r) = 
8(oo)—6(—oo) = 2n. It is well known that the instanton (anti-instanton) is almost universally 
constant except for a very small region of time variation. In the QPS process the amplitude 
of the superconducting order parameter is zero only in a very small region of space as a 
function of time and the phase changes by ±2n. So our system reduces to a neutral system 
consisting of equal number of instanton and anti-instanton. One can find the expression for 
0(u) after the Fourier transform applied to both the sides of Eq. 11 which yields 0(u>) = 
z J2i eih(iu)e tuJT1 . Now we substitute this expression for 9{uj) in the second term of Eq. 9 



and finally we get this term as J2ij F{ T % ~ T j) e i e j, where Ffa — Tj) = ^ J2 m y^-re"^ -7 "^ 
~ ln{ji — Tj). We obtain this expression for very small values of co ( — > 0). So F{ji — Tj) 
effectively represents the Coulomb interaction between the instanton and anti-instanton. 
This term is the main source of dissipation of SQUID array system. Following the standard 
prescription of in^tirne^ath integra. formal, we can write the partition function 



of the system as [8j, |10| 
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Z = E-T^E / dT n / dTn-1- 

^ n\ ^Jo Jo 

r dT ie - F{n - T ^ . (12) 

Jo 

We would like to express the partition function in terms of integration over auxiliary field, 
q(r). After some extensive analytical calculations, we get 

Z = J Z} g (r)e ( "S»"^ 5(!W " )9( " iu " ))+(22 'C' lrc ' )S ' !(T)) . (13) 

Thus by comparing the first term of the action of Eq. (8) and the first term of exponential 
of Eq. (13), we conclude that the dissipative strength a and the Luttinger liquid parameter 
of the system are related by the relation, K = 4a. 



3. QUANTUM FIELD THEORETICAL STUDY OF MODEL HAMILTONIAN OF 
THE SYSTEM AND EXPLICIT DERIVATION OF DISSIPATIVE STRENGTH 

In the previous section and also from our previous study [10(], we have shown explicitly 
that the mesoscopic SQUIDs array is equivalent to the array of superconducting quantum 
dots (SQD) array with modulated Josephson coupling. We first write the model Hamiltonian 
of SQD with nearest neighbor (NN) Josephson coupling Hamiltonian {Hj{) and also with 
the presence of the on-site (Heco) and NN charging energy (Heci) between SQD as, 

H = Hji + H E co + H E ci- (14) 

We would like to recast our model Hamiltonian in the magnetic flux induced Coulomb 
blocked regime (Eco » Ej). In this regime, one can recast the Hamiltonian in spin 
operators. It is also observed from the experiments that the quantum critical point exists 
for larger values of the magnetic field, when the magnetic field induced Coulomb blockade 
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FIG. 2: Quantum phase diagram {Ej vs. V g ) of the SQD array. We have depicted the different 
phases of the model Hamiltonians by: A. Mott insulating phase, B. Charge-density wave (CDW); 
C. First kind of Repulsive Luttinger liquid (RL1); D. Second kind of Repulsive Luttinger liquid 
(RL2); E. Superconductivity. Q and P2 are the multi-critical points (please see the text). K is 1 at 
the phase boundary between Luttinger liquid (RL2) and superconductivity and also at the phase 
boundary between Mott phase and superconductivity. K = 1/4 at the phase boundary between 
RL1 and CDW state. K = 2 and 1/2 at the PI and Q point respectively. Q and PI are the special 
points, Q is the charge degeneracy point and PI is the particle-hole symmetric point. 



phase is more prominent than the Ej induced SC phase. Thus our theoretical model is 
consistent with the experimental findings. During this mapping process we follow Ref. (Q- 
Q), so that. 

Hjt = -2Eji'£{SfS i+1 - + h.c) 



E n = EAcos 



7T$ % 
*0 ' 



H 



E, 



ECO 



CO 



E(s t z -h) 2 . 



are the Hamiltonian Heco accounts for the influence of gate voltage (eiV ~ V^), where eN 
is the average dot charge induced by the gate voltage. When the ratio — > 0, the SQD 
array is in the insulating state having a gap of the width ~ Eco, since it costs an energy 
~ Eco to change the number of pairs at any dot. The exceptions are the discrete points at 
iV = (2n + 1), where a dot with charge 2ne and 2(n + l)e has the same energy because the 
gate charge compensates the charges of extra Cooper pair in the dot. On this degeneracy 




FIG. 3: Shift of the Josephson coupling (AEj) due to the magnetic flux induced dissipation effect 
with magnetic flux (measured with respect to magnetic flux quantum, 4>o = ||). The red line 
is for the shifting of the Josephson coupling for the phase boundary between the RL2 and the 
superconducting phase. Blue line is for the shift of the phase boundary between the CDW and 
RL1 phase. Inset shows the shift of the Josephson coupling for the particle-hole symmetric point 
(green line) and the Charge degeneracy point (yellow line) in the figure. 

point, a small amount of Josephson coupling leads the system to the superconducting state. 
Here h = N ~ 2 ^ 1 allows the tuning of the system around the degeneracy point by means 
of gate voltage. H EC i = ^E Z \ S z Si + \ Z . At the Coulomb blocked regime, the higher 
order expansion leads to the virtual state with energies exceeding Eqq- In this second order 
process, the effective Hamiltonian reduces to the subspace of charges and 2, and takes the 



form 



He = -^Y.Si Z S^ z - ^US^Sr + h.c). (15) 



With this corrections H EC i become 



3E 



2 



H EC1 ~ {AE Z1 - J2S z S i+1 



4K 



CO 



In this analytical expression, we only consider the nearest -neighbor contribution of the inter- 
action. There is no evidence of next-nearest-neighbour interaction for mesoscopic SQUID ar- 



ray system 13j. One can express spin chain systems to as spinless fermions systems through 
the application of Jordan- Wigner transformation. In Jordan- Wigner transformation the re- 
lation between the spin and the electron creation and annihilation operators are S z (x) = 
ijj^(x)ip(x) — 1/2, S~(x) = ip(x) exp[i7r X^=ioo %]> where S + = (S~) + , n(x) = tfj" f (x)ip(x) 
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FIG. 4: Dissipative quantum phase diagram (Ej vs. V g ) of the SQD array. We have depicted the 
different phases of the model Hamiltonians by: A. Mott insulating phase, B. Charge-density wave 
(CDW); C. First kind of Repulsive Luttinger liquid (RL1); D. Second kind of Repulsive Luttinger 
liquid (RL2); E. Superconductivity. Q and P2 are the multi-critical points (please see the text). 
K is 1 at the phase boundary between Luttinger liquid (RL2) and superconductivity and also at 
the phase boundary between Mott phase and superconductivity. K = 1/4 at the phase boundary 
between RL1 and CDW state. K = 2 and 1/2 at the PI and Q point respectively. Q and PI 
are the special points, Q is the charge degeneracy point and PI is the particle-hole symmetric 
point. We notice that due to the dissipation, phase boundaries and special points are shifted w.r.t 
Josephson coupling, please see the text for detail explanation. 

is the fermion number at the site x. We have transformed all Hamiltonians in spinless 
fermions as follows: H n = - 2 E n E^iVj+i" + h.c), H E co = 2hE C o EiCv^i - V 2 )- 
H ECl * (4£ Z i - 3 ~M) iMHi ~ ~ V 2 )- 



In order to study the continuum field theory of these Hamiltonians, we recast the spinless 



30). 



fermions operators in terms of field operators by a relation 

^(ar) = [ e ikFX ^ R {x) + e~ lkpX ^ L {x)\ (16) 
where ipn(x) and iPl(x) describe the second-quantized fields of right- and the left-moving 

11 



fermions respectively. We would like to express the fermionic fields in terms of bosonic field 
by the relation 

^ r ( x ) = JL= e -i WW - W), (17) 
y2ixa 

where r denotes the chirality of the fermionic fields, right (1) or left movers (-1). The 
operators U r preserve the anti-commutivity of fermionic fields. field corresponds to the 
quantum fluctuations (bosonic) of spin and 9 is the dual field of 0. They are related by the 
relations 4>r = 9 — <fi and <p L = 9 + 0. The Hamiltonians without (Hi) and with 
co-tunneling (H 2 ) effect are the following 



Hi = H + ^ / dx : cos(A\fK<j)(x)) : 

(2ira) J 

J (d x cf>(x)) dx (18) 



(4E~ - 



H 2 = H + §^ f dx : cos^Vk^x)) : 

(2na) J 

+— / (dA(x)) dx (19) 

Where, H is the non-interacting part of the Hamiltonian. The Luttinger liquid param- 
eters of the Hamiltonian Hi and H 2 are Ki and K 2 are respectively. 

7T + 2sin~ 1 Ai v 1 

K 2 = 7T- (21) 

Ti + 2sin~ 1 A 2 v ; 

Ai = A 2 = — §§^- We calculate the dissipation strength by calculating K 

for both cases and then we use the relation K = 4a. Therefore the dissipative strength in 
absence (ai) and presence (a 2 ) of co-tunneling effect are 

«i = nr~ (22) 

012 = L + 2 L-ia 2 (23) 

This is the first analytical derivation of flux induced dissipation strength in terms of the 
interactions of the system. In the derivation of K 2 and a 2 , we only consider the nearest- 
neighbour hopping consideration when we consider the effect co-tunneling. We consider these 
two processes to emphasis the importance of co-tunneling effect for this system. Before we 
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proceed further for the analysis of the Hamiltonian , H\ and if 2, we want to explain in 
detail for the different values of K at the different points (like PI, Q and P2) of phase 
diagram and also at the phase boundaries. The physical analysis of the phases and clear 
distinction of the phase boundaries will depend on the values of the K. Point Q is the 
charge degeneracy point for the low Cooper-pair density (n, = 1/2). The value of K will be 
evaluated from the relevance of sine-Gordon terms. At the point Q, system is in the second 
order commensurability, sine-Gordon coupling terms ,Eq. 18 and Eq. 19, will become 
relevant for K = 1/2, which is depicted in the Fig.l. Point PI is the particle hole symmetric 
point, i.e., there is one particle in each site. At this point system is in the first order 
commensurability, i.e., the sine-Gordon coupling term is cos(2\/K(f)(x)). So this term will 
become relevant for K=2, which is depicted in the Fig.l. We are understanding from Eq. 18 
and Eq. 19 that the applied gate voltage acts as a chemical potential. So the proper tuning 
of gate voltage will drive the system from the insulating state to the other quantum phases 
of the system. We are now interested in finding the value of K at the phase boundaries, 
hence analysis is the following: We follow the Luther-Emery [29] trick during the analysis. 
One can write the sine-Gordon Hamiltonian for arbitrary commensurability as 

H 3 = H + \J dx cos(2nVK(f>(x)), (24) 

where n is the commensurability and A is the coupling strength. H is 
the free part of the Hamiltonian. We know that for the spinless fermions, 
V'-rVl + ipL ipR — 2™^ / dxcos{2\fK(j){x)) , which is similar to the analytical expres- 
sion of sine-Gordon coupling term but with the wrong coefficient inside the cosine. One can 
set 4>{x) = 2\/~K(j)(x) then the Eq. 24 become 

H A = H + \J dxcos(2(p(x)). (25) 

K and K are related by the relation, K = 

At the phase boundary, K — 1 that implies K = 1/n 2 . So for the first and second order 
commensurability the value of K at the phase boundary are 1 and 1/4 respectively which 
is depicted in Fig. 1. The point to be noticed that if we start from an initial model with 
K = 1/n 2 , i.e., in general a strongly interacting model, the resulting spin-less fermions 
model corresponds in the boson language to K = 1 which means that it is non-interacting. 
For this particular value of K the spin- less fermions whose bosonized form is Eq. 25 are just 
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ree particle with backscattering. This special value of K is known as the Luther-Emery 



29j | line, the importance of the Luther-Emery solution is to provide a solution for the 



massive phase on the whole line K = 1/n 2 for arbitrary A. 

Here we do the analysis for Hamiltonian Hi. In the limit A = Ai, for Ej < 2Ezi 
and relatively small field, the anti-ferromagnetic Ising interaction dominate the physics of 
anisotropic Heisenberg chain. When the field is large, i.e., the applied gate voltage is large, 
the chain state is in the ferromagnetic state. In the language of interacting bosons, The 
Neel phase is the commensurate charge density wave phase with period 2, i.e., there is 
only one boson in every two sites. In Fig. 1 this phase region is described by region B. 
The ferromagnetic state is the Mott insulating state, this is the phase A of our quantum 
phase diagram (Fig. 1). Hi is the Heisenberg XXZ model Hamiltonian in a magnetic 
field. The emergence of two Luttinger liquid phases for the following reasons: RL1 and 
RL2 respectively occur due to commensurate-incommensurate transition and the criticality 
of Heisenberg XY model. For the intermediate values of the field, system is either in the 
first kind of repulsive Luttinger liquid (RL1) phase for K < 1/2 or in the second kind 
of repulsive Luttinger liquid (RL2) phase for K > 1/2. The physical significance of RL1 
phase is that the coupling term is relevant but the applied magnetic field, i.e., the applied 
gate voltage on the dot, breaks the gapped phase whereas in the RL2 phase non of the 
coupling term is relevant due to the larger values of K (> 1/2). The phase regions are 
described by C and D respectively for the RL1 and RL2 in Fig. 2. These phase regions in 
Fig. 2, are all most same for two different limits, A = Ai and A = A2. So we predict the 
existence of two RL from two different sources. In previous studies this clarification was 
absent and they had reported only one RL [8j. The value of K = 1 at the phase boundary 
between the MI and SC phase and also between the RL2 and SC phase. From the analysis 
of K at the phase boundary we obtain Ez\ = 0, according to our theory and also from the 
experimental findings this condition is unphysical. So the interaction space of Hamiltonian 
Hi is not sufficient to produce the whole phase diagram of Fig. 2. It indicates that we 
shall have to consider more extended interaction space to get the correct phase diagram. 
If we consider the co-tunneling effect in this Hamiltonian system, i.e., the Hamiltonian 
H 2 . The phase boundary analysis at the MI and SC phase and also for the RL2 and SC 



phase implies that we get the condition Eji = ^IMLziMco^ which is consistent physically. 
Now we discuss the effect of magnetic field induced dissipation on the quantum phase 
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diagram of superconducting quantum dot lattice. We have already proven in the previous 
section the analytical relation between the LL parameter (K) and the dissipation strength 
in presence of magnetic flux. Here we use the analytical expression for K to study the 
effect of magnetic flux on the quantum phases and quantum phase boundaries. Therefore 
the modified quantum phase diagram in presence of magnetic flux include the magnetic 
flux induce dissipation effect. In the previous paragraph we emphasis the importance of 
co-tunneling effect. Therefore we consider the Eq.(21) and Eq.(23) when we consider the 
shift of the Josephson couplings due to the presence of magnetic flux induced dissipation. It 
is very clear from Eq. 21 and Eq. 23 that the analytical relation of the Josephson coupling 
with the system interactions parameters at different phase boundaries and different points 



are more complicated than the Ambegaokar and Baratoff relation j27|. The analytical 
expression for the Josephson couplings at the different phase boundaries and special points 
are the following: 

Particle- Hole symmetric point (K = 2), 



Ej = 0.9428£coG/l + + 1) (26) 
V &co 



Charge-Degeneracy point (K = 1/2), 



Ej = 0.75^/1 + ^-1) (27) 



Phase boundary between the RL2 and SC phase (K = 1), 



E z %Ecqs 



Ej = 4(W-^-^) (28) 



Phase boundary between CDW state RL1 phase (K = 1/4), 



AEj = 0.75E C0 (Jl + ^ £i + l)- (29) 
V ^co 

Our main intension is to study the effect of magnetic flux on the particle hole symmetric 
point (K = 2 ), charge degeneracy point (K = 1), multicritical point and also for the phase 
boundaries between the different quantum phases. 

The analytical expressions for the shift of Josephson couplings (AEj) due to the presence 
of magnetic flux are the following: 
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Particle- Hole symmetric point (a = 1/2), 



AEj = 0.9428^00(4/1 + ^ + l)(l/\cosA\ - 1) (30) 
V &co <Po 



Charge- Degeneracy point (a = 1/8), 



AEj = 0.75E C0 (Jl+ 3 -§^-l)(l/\cosA\-l) (31) 
Phase boundary between the RL2 and SC phase (a = 1/4), 



AEj = 4(^^)(l/|co S (|)| - 1) (32) 



Phase boundary between CDW state RL1 phase (a = 1/16), 



AEj = 0.75^0(1/1 + ^ + l)(l/\cosA\ - 1) (33) 
V &C0 <Po 

In Fig. 3, we present our results of deviation of Josephson coupling as a function of 
magnetic flux. We observe that this deviation is maximum at the half-integers values of 
magnetic flux quantum. It reveals from our study that the effect of magnetic flux induce 
dissipation for the half-integer values of magnetic flux quantum is slightly prominant for 
the special points compare to the qunatum phase boundaries. Inset shows the shift of 
Josephson coupling for the particle-hole symmetric point and the charge degeneracy point. 
It is also clear from the inset that the magnetic flux induce dissipation is more prominant 
for the particle-hole symmetric point compare to the charge degeneracy point. 
It is very clear from the analytical expression from Eq. 30 to Eq. 33 and also from the 
Fig. 3 from our study that there is no appreciable changes in the quantum phase diagram 
for small magnetic flux in the system. As the applied magnetic fluxes changes from O.40 o 
to 0.600) quantum phase diagram shows some appreciable change and it is robust for the 
half-integer magnetic flux quantum. 

In Fig. 4, we present the magnetic flux induce dissipative quantum phase diagram of our 
system. This schematic phase diagram is for values of magnetic flux which appreciably effect 
the phase boundaries and special points as we have discussed in the previous paragraph. 
We present the phase boundaries and special points in terms of the dissipative strength of 
the system. We observe from our study that magnetic flux induce dissipation favour the 
insulating phase and the gapless LL phase over the superconducting phase of the system 
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which is consistent with the experimental findings [13]. 

SUMMARY AND CONCLUSIONS 

We have studied the quantum phase diagram of mesoscopic SQUID array in absence 
and presence of magnetic flux. The magnetic flux induced dissipation modified quantum 
phase diagram of our system. We have derived an analytical relation between the Luttinger 
liquid parameter and dissipation strength. We have also noticed that magnetic flux induced 
dissipation effect is not same for all values of magnetic flux quantum. We have also observed 
that the magnetic flux induce dissipation favours the insulating phase of the system over 
the Luttinger liquid and superconducting phase of the system. 

Acknowledgement: The author would like to acknowledge Dr. N. Sundaram for reading the 
manuscript critically. 
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